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We consider thermodynamics of the excluded volume particles at finite temperature and chemical 
potential, in the low density approximation. We assume Boltzmann statistics and study the influence 
of the excluded volume on an ideal gas thermodynamics at the same temperature, pressure and 
numbers of particles. We show, that considering the change of the free enthalpy due to the excluded 
volume, and using the Maxwell identities, one can derive relevant thermodynamic functions and 
parameters of multi-component gases. The derivation is quite general as particles may have different 
sizes and shapes which can also depend on their momenta. Besides it’s simplicity and generality, 
our approach has the advantage of eliminating the transcendental equations occurring in earlier 
studies. A representative example of the excluded volume thermodynamics is the single-component 
gas of hard spheres. For this case, using the viral expansion, the validity limits of the low-density 
approximation are also discussed. 


I. INTRODUCTION 

A detailed analysis of secondary hadrons produced in 
heavy ion collisions (HIC) in a very broad energy range 
from SIS to LHC has shown that they are of thermal ori¬ 
gin M. Particle yields are excellently described as the 
Hadron Resonance Gas (HRG), an uncorrelated gas com¬ 
posed of all known hadrons and resonances constrained 
by the conservation laws [l|. The HRG is also very suc¬ 
cessful to quantify the equation of state (EoS) of the 
hadronic phase of QCD at finite temperature and den¬ 
sity, as was shown recently by comparing the EoS ob¬ 
tained in lattice gauge theory with the predictions of the 
HRG model [^-Q- These results indicate that HRG is a 
remarkably good approximation to QGD thermodynam¬ 
ics in the confined phase, and that HIG experiments are 
indeed probing the thermal QGD medium. 

The statistical partition function for HRG is con¬ 
structed for a mixture of ideal gases, nevertheless, ac¬ 
cording to the arguments by Dashen, Ma and Bernstein 
3 and Hagedorn [T^ , it includes attractive interactions 
through resonance contributions. 

Different extensions of the HRG model have been pro¬ 
posed, including not only the attractive, but also the re¬ 
pulsive interactions between hadrons. The later are ex¬ 
perimentally known to be important at short distances. 
Glearly, the effect of both interactions on the thermody¬ 
namics can be introduced in a model independent way, 
using the S-matrix approach Q. In the low density ap¬ 
proximation the resulting expression for the partition 
sum of the hadron gas reduces to the Beth-Uhlenbeck 
form [Tlj| . where the interaction contribution is linked 
to the experimental phase shifts dl. However, due to 
far incomplete data on phase shifts, this approach can 
be systematically used only in selected sectors of hadron 
interactions [l3 - [l^ . 

A transparent phenomenological way to account for 
repulsions in a hadron gas is to abandon the point-like 


nature of particles and introduce for them finite sizes. In 
this way the repulsive interactions in the hadron gas are 
implemented via the excluded volumes. Exact formulae 
for the thermodynamic functions of the excluded volume 
gas are not known, but low density approximations have 
been developed (isl - l^ . The hard core hadron resonance 
gas models were also successfully applied in heavy ion 
phenomenology [iHsll and in the interpretation of recent 
lattice data |26l - l^ . 

In the following, we focus on the implementation of 
the repulsive interactions in a hadro n g as via excluded 
volumes, as introduced in Refs. [^[^. We propose a 
new method, which is based on the Maxwell identities, 
to calculate the thermodynamic functions and parame¬ 
ters of the low density excluded volume gas, subject to 
Boltzmann statistics. They are expressed in terms of the 
excluded volumes for pairs of particles, and the thermo¬ 
dynamic functions and parameters of the ideal gas at the 
same temperature, pressure and number of particles of 
each kind. The derivation is quite general, as particles 
may have different sizes and shapes, which may also de¬ 
pend on their momenta. Besides its simplicity and gen¬ 
erality, our approach has the advantage of eliminating 
the transcendental equations occurring in earlier studies. 
The results from Refs. [20l-l^ follow as special cases 
from our analysis, except for the partial pressures in an 
r-component gas of hard spheres presented in Ref. , 
where some additional assumptions, or the higher order 
corrections, are needed. 

A system of hard spheres with all particle volumes be¬ 
ing equal, is a special example of the excluded volume 
gas. For this case we use the known results concerning 
the virial expansion, to find the validity limits of the low 
density approximation to the thermodynamics of the ex¬ 
cluded volume gas. 
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II. THERMODYNAMICS OF THE EXCLUDED 
VOLUME GAS 

When modeling thermodynamics of the extended par¬ 
ticles, a dimensional scale is the size of the excluded vol¬ 
umes in coordinate space, which is the input in the phe¬ 
nomenological approach, whereas a basic variable is the 
configuration integral. 


For the case of N identical particles with the unique 
two-particle excluded volume 2v(T), one finds,that in 
the large V-limit 

f V N \ 

r log C = NT i^log - - -v{T) + 1 j , (7) 

and 


C = 


j3Ar 


xe 




( 1 ) 


pV 

Wf 


V dlogC 
N dV 


N,T = l+v{T)n, 


( 8 ) 


where the summation is taken over all pairs of particles. 
A pair-wise potential in the exponent, is introduced such 
that, Uij{xi — Xj) = oo when particles i and j overlap, 
and Uij{xi — Xj) = 0 otherwise. The occurrence of par¬ 
ticle overlap may depend on more variables than written 
explicitly in Eq. m, e.g. it can depend on the orienta¬ 
tions of the volumes associated with the particles. Then 
it is understood, that in Eq. o, the integrations over 
these additional variables are also included. 

The definition given in Eq. o is very general, as it is 
applicable for any shapes of the particles, which are not 
necessarily the same for all of them. The only require¬ 
ment is, that the pair-wise potential must be define when 
two given particles overlap and when they do not. 

Introducing the Mayer functions, 

( 2 ) 

in Eq. o, expanding the integrand of C in powers of 
fij up to leading order, and performing the integrations, 
one finds, that 


where 



( 3 ) 

( 4 ) 


are the excluded volumes for pairs of particles, which in 
general, may depend on temperature. 

The evaluation of the integrals over Mayer functions in 
Eq. dH) can be elementary, like e.g. for the gas of hard 
spheres, or may be more advanced, like for the gas of 
Lorentz contracted hard spheres [l^. In the following, we 
will not introduce any particular form of fij , and assume, 
that the integration in Eq. (U) has been done. 

For further applications we consider the expression, 
TlogC. From Eq. ([3]), one finds, that to leading or¬ 
der in Vij, 


T 


TlogC = NT\ogV - AT)- 


V 

{iJ} 

This yields, in particular, the equation of state 

N 


( 5 ) 


pV VdlogC, 
NT~ N dV 


k.T = 1 + ^ (6) 


{ij} 


where n=N/V, is the particle density. Note, that when 
the particles are not identical, but all Vij{T) are equal, 
then only the combinatorial factor in the configuration 
integral changes, so that the equation of state remains 
unchanged. 

In some approaches, the excluded volume depends on 
the momenta of the particles (cf. e.g. [l^). Then, the 
Mayer function for particles i and j depends on their 
momenta and it is necessary to average it over the Boltz¬ 
mann distribution. As a result, the coefficients Vij be¬ 
come temperature dependent. Having made this change, 
the terms of first order in the excluded volumes are cor¬ 
rectly reproduced. However, since the average of the 
square of the excluded volume is not equal to the square 
of its average, this simple method works only up to first 
order. This is enough for the low density approxima¬ 
tion, however not for the more precise calculations needed 
when discussing its validity limits. 

Let us consider an r-component gas and denote by Ni 
the number of particles of species i, with the total number 
of particles YA.=i ~ Then, including the combina¬ 
torial factor, 01=1 Tfif in the configuration integral o, 
one finds, that in the large A^i-limit for all i = 1,... r, 

TlogC = NTlogV - (^) (9) 

- T^iV,(log7V,-l), 


and 


= i + (10) 

^,3 

The above result agrees with that obtained in Ref. [^ . 
However, contrary to the earlier approaches, there is no 
need to solve any transcendental equations to describe 
relevant thermodynamic functions. The reason is that 
we use a different set of parameters as independent vari¬ 
ables. For the description of a state of the system any 
complete set of parameters is equally acceptable. E.g. 


For N identical particles, one needs to introduce the combina¬ 
torial factor, 1/A!, into the configurational integral fTIl . 
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for the single-component gas, one can use the parame¬ 
ters {T,p,N) or {T,V,p). The fact, that in some pro¬ 
cess N changes at constant /i, or the other way round, 
is irrelevant here. Indeed, the pressure, temperature and 
chemical potential are related by the following equation 

p = (11) 

In Ref. [ 2 ^ this equation is used to find p{T,fj,), which 
requires the solution of a transcendental equation. The 
calculation of p{p,T), from Eq. dTlT) is, on the other 
hand, elementary. 


III. MAXWELL IDENTITIES AND THEIR 
IMPLICATIONS 


Elementary phenomenological thermodynamics will 
now be used to find the thermodynamic functions and 
parameters of the excluded volume gas at low density, 
in terms of the excluded volumes Vij(T), and the pa¬ 
rameters and thermodynamic functions of an ideal gas 
at the same temperature, pressure and numbers of par¬ 
ticles. The parameters and thermodynamic functions of 
the ideal gas will be denoted with the superscript, id. 

Let us consider first a single-component gas of particles 
with the excluded volume for a pair of particles, 2v(T). 
In this case, the parameter v{T) may be interpreted as 
the single-particle excluded volume. Now we introduce 
the excluded volume 2Xv(T). Then for A = 0, we have an 
ideal gas, and for A = 1, the gas we are interested in. The 
point is to find, how the thermodynamic functions and 
parameters change, when A increases from zero, where 
the familiar relations for an ideal gas hold, to one. The 
changes are calculated at constant temperature, pressure 
and number of particles. 

We consider the differential of the free enthalpy 

dG =—SdT + Vdp + pdN + pNv{T)dX, (12) 


where the last term on the right-hand side is the work 
done on the system when the excluded volume of each 
particle grows. This work is done against the pressure 
in the surrounding of the particle, which is equal to the 
external pressure p. The low density approximation is 
implied by the form of this term. In the exact theory, the 
excluded volumes can overlap and, therefore, the increase 
of the total excluded volume is slower than linear in N. 
For a low density gas, however, these overlaps can be 
neglected and formula (TT^ holds. 

Since free enthalpy is a well-defined function of state, 
the following Maxwell identities are valid 




\N,T,p = 


\N,T,p = 


\N,T,p = 


dpNv{T) 

dpNv{T) 

dp 

dpNv{T) 

Wf 


T.p,v=pv{T), (13) 

T,N,v = Nv{T), (14) 

dv{T) 


N,p,v — pN - 


dT 


(15) 


Integrating the above equations at constant temperature, 
pressure and number of particles from A = 0 to A = 1, 
one gets 


p{T,p) 

= p^‘^iT,p)+pv{T) 

(16) 

V{T,p,N) 

= V^‘^{T,p,N) + Nv{T), 

(17) 

S{T,p,N) 

= S^‘^iT,p,N)-Np^^. 

(18) 

For the energy, this implies 


EiT,p,N) 

= E^d^T,p,N)-TNp^^. 

(19) 


Equivalent expressions for thermodynamic functions 
have been obtained in Ref. [l^, with the substitution, 
NT = pV, in the correction term in Eq. CSl), which is 
legitimate at first order in v(T). 

For a Bose gas the well-known condition, p,^‘^{T,p) < 
m, where m is the mass of a particle, implies that 

p{T,p) < m + pv. (20) 


This upper bound is equal to the minimum free enthalpy 
needed to introduce one more particle into the system. 
It consists of the rest energy of the particle m and of 
the work pviT), at constant temperature and pressure, 
necessary to create the single particle excluded volume. 
This observation implies, that when v(T) decreases, the 
abundance of particles is enhanced. 

Using the expressions for the parameters V, S, E of the 
excluded volume gas in terms of the parameters of the 
ideal gas from Eqs. [17,18,19], one finds the relations 
between the volume densities, 


n 


s 


e 


^id 


1 -I- 

id _ id dv{T) 

^ P dT 

1 -I- vrd‘^ 

Td - 


1 -I- vri 


id 


( 21 ) 

( 22 ) 

(23) 


of particle number, entropy and energy, respectively. The 
arguments (T, p, N ), have been skipped in these equa¬ 
tions. 

The above relations, which link the thermodynamics of 
the point-like and the excluded volume particle gases, are 
consistent with previous studies in Ref. [ 2 ^ , and are also 
very successful in the applications to thermal description 
of particle production yields in HIC [H-Q . 

The method introduced above, for the single¬ 
component gas, can be generalized to multi-component 
system. For an r-component gas the differential of the 
free enthalpy reads, 


r 

dG = -SdT + VdpT'^p^dNi (24) 
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where we have introduced, pi = -jfp, the relation which 
is valid for ideal gases, and is also legitimate, in the cor¬ 
rection term due to excluded volumes. 

The last term in Eq. (IMl) . which describes the work 
done when the excluded volumes increase, can be inter¬ 
preted as the sum of works done around all particles. 
The work done around a particle of type i is a sum of r 
terms: the work needed to free the volume Vi^i{T)d\ from 
particles of type one against the partial pressure pi , the 
work needed to free the volume Vi^ 2 {T)dX from particles 
of type two against the partial pressure p 2 , and so on, up 
to particles of type r. Thus, for the multi-component gas 
the one-particle excluded volume can be defined, but it 
does not have a simple geometric meaning. 


and numbers of particles, one finds 


d 

V 

s 


id 1 XliNj 

* j 


yl. 


3 ^ 

1 iV^ 
2 ^^ N 




Vi- 




1 N^Nj dv-i^j 

~/V 


(28) 

(29) 

(30) 


where the arguments in each case are {T,p, Ni,..., Nr). 
Following, the discussion of the single-component gas, 
the generalization of Eqs. (1^ . (1^ and (1^ . to the 
r—compannet gas, is rather transparent. 


Moreover, there is some ambiguity in dividing the 
work associated with the increase of the excluded vol¬ 
umes among the surroundings of the particles. Clearly, 
the term pi 7 V 2 Wi 2 (iA=p 2 .^i'Ci 2 dA, contributes to the work 
done in the surroundings of the particles of types 1 and 
2. Following the interpretation presented above, for each 
particle of type 1 this work is ^Vi 2 dX, and for each par - 
ticle of type 2, it is ^vi 2 dX. According to Ref. j^ . 
however, for r = 2 and for the second kind of particles 
being point-like, the contribution of this term to the work 
done around a particle of type 1 is zero, and around a 
particle of type 2 it is piVi 2 {T)dX. Both these interpre¬ 
tations are consistent with our leading order expression 
for the differential of free enthalpy. 

For an r—component gas, the Maxwell identities are 
generalized as 


dp.i I 






dN, 


N. 




iT.p.A 


(25) 


2"^ 7V2 


^1,3^ 


a(li'E.j^*'.,i(r)) 




dp 


\N,T,X 


N,N. 


2^ N 




(26) 




dT 

1 NiNj dvij{T) 

o 2^ 


\N,p,X 


2 ^ N 


dT 


(27) 


where the subscript N indicates that all the multiplicities 
Ni ,..., Nr are fixed. Integrating the above equations 
over A from zero to one, at fixed temperature, pressure 


IV. VALIDITY OF THE LOW DENSITY 
APPROXIMATION 

In the previous sections, we have derived the thermo¬ 
dynamic observables of a gas of extended particles, under 
the low density approximation. In order to determine re¬ 
liably the validity limits of this approximation, it is nec¬ 
essary to compare our results with the exact solution, or 
with its good approximation. In the excluded volume gas 
problem, such comparison is possible for the gas of hard 
spheres of the same radius i?, i.e. with the volume 

The corresponding single-particle excluded volume is 
then, V = Avq. For R = 0, the gas reduces to the ideal 
gas, and the low density approximation yields the ex¬ 
act result. With increasing R the approximation dete¬ 
riorates. The task is to find, up to what values of R, 
the low density approximation is reliable, within a given 
error margin. 

The equation of state of the gas of hard spheres can be 
written in the form of a virial expansion, 

OO 

^=EfciZ*"^ z = nv. (32) 

The virial coefficients bi are known explicitly [H, [s^ for 
i = I,..., 12, and they are quoted in Table I. 


bi 

1 

64 

0.286949 

by 

0.0130235 

^10 

0.0004035 

b2 

1 

^5 

0.11025 

bs 

0.004183 

bn 

0.000123 

bs 

0.625 

be 

0.0388819 

bg 

0.001309 

bi2 

0.000037 


Table I. Virial coefficients for the gas of hard spheres 
from Refs. [33.13^. 

Comparing Eqs. (ED and (l32l) . it is clear, that the low 
density approximation is obtained when only the first two 
terms of the virial expansion are taken into account. The 
convergence of the virial expansion deteriorates when the 
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FIG. 1: Left-hand figure: The virial series from Eq. dMl) truncated at m = 2,..., 12. The term m = 2 corresponds to the 
lowest line. Right-hand figure: The limiting values zum of z for which the accuracy of Sm{z) from Eq. (I5S1) with m=2,..,10, 
relative to Si 2 {z), is 2% (diamonds) and 10% (dotes). 


parameter ^ increases. In order to find how good the 
convergence is, we must define the relevant range of z. 

Obviously, the total volume of the spheres Nvq must 
be smaller than the total volume of the gas. This implies, 
that 

z<A. (33) 

Moreover, it is not possible to pack hard spheres into a 
volume without leaving free spaces. Consequently for the 
closest packing 

An equivalent formula was conjectured by Kepler in 1611, 
then Gauss in 1831 demonstrated, that this is the clos¬ 
est packing possible, if the spheres form a lattice. The 
general proof is difficult, however can be handled by us¬ 
ing computers. Its completion was announced by T.C. 
Hales in 2014. From this limit, one gets 

2 < 2.96. (35) 

Finally, it is plausible, that the hadrons form a fluid 
rather than a condensed phase. The condition for the 
freezing of the gas of hard spheres is, nuo = 0.494 [s^, 
which yields the limit, 

z < 1.976. (36) 

Above this limit, the virial expansion for the gaseous 
phase of hard spheres system diverges. Consequently, one 
concludes, that the relevant region of z, is as follows 

0 < z < 2. (37) 

We introduce the partial sum, 

m 

Sm{z) = Y,hz^-\ (38) 

i=l 


for m > 2. Knowing the values of virial coefficients, one 
can quantify, how good is the low density approximation 
for different values of 0 < z < 2. 

In Fig. [TJleft, we show the sum Sm{z) for 0 < z < 2, 
and for different m = 2,..., 12. The low density approx¬ 
imation is seen in this figure to describe the equation of 
state within less than few percent, for z < 0.5. With 
increasing z up to unity, the series is known to converge 
and the deviations are up to 30%, whereas for z > 1 a 
large corrections to the low density approximation are to 
be expected. For z close to two, the virial expansion in 
Eq. (1^ diverges. 

In Fig. [T]-right, we show for m = 2,... 9 the limiting 
values of z for which the accuracy of approximation of 
Si 2 iz) series by Sm{z) is two and ten percent, respec¬ 
tively. Since only the values of z < 2 are of interest, the 
limits significantly exceeding z = 2, are not shown in 
this figure. 

The Hadron Resonance Gas (HRG) formulated with 
excluded volume from Eq. m, was applied to describe 
particle production in heavy ion collisions in a very broad 
energy range from SIS up to LHG [I| . All particle yields 
in heavy ion collisions were successfully quantified by the 
model along a common freezeout line in the tempera¬ 
ture and baryon-chemical potential plane [ms. This 
finding allows to quantify the conditions for the appli¬ 
cability of the low density approximation in the HRG 
of excluded volume particles, for the phenomenologically 
relevant system, along the freezeout line. 

In applications, the total density of particles at chem¬ 
ical freezeout is temperature and chemical potential de¬ 
pendent, and consequently, changes with the collision en¬ 
ergy ^/s. Thus, z in Eq. (l32|) is a function of ^/s, and of 
the chosen value of R. Let us introduce the function 

z{'/s,R) = (39) 

where R is expressed in fermis. 

In Fig. [5]-left, we show the energy dependence of 
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FIG. 2: Left-hand figure: Function ^(-^5, 0.3) from Eq. (|39ll along the chemical freezout line extracted from particle yields 
data taken in central nucleus-nucleus collisions at different collision energies P, |M|- Right-hand figure: The limiting values of 
R from Eq. (I39II at n = 0.375fm~®, for which the accuracy of the truncated series Sm{z) for m=2,..,10, relative to Si 2 {z), is 
2% (diamonds) and 10% (dots). 


z(y/s, 0.3) for the HRG, which includes all particles and 
resonances, listed by the Particle data Group. For sim¬ 
plicity, the radius R in Eq. (15^ was assumed to be com¬ 
mon for all particles and antiparticles with the value, 
R ~ 0.3 fm. 

The function 2 :(y^, 0.3) is small at low energies bel¬ 
low AGS, increases monotonically with increasing energy, 
and at high energies beyond SPS, it saturates at the 
value slightly below 0.17. In our further estimates we 
will use the approximation z(-y/s, 0.3) = 0.17, or equiva¬ 
lently, n = 0.375 fm“^. 

The upper bound of z, established in Eq. (I36p implies, 
in general, the restriction on the radius of the hard core 
repulsion, 

R < 0.680 fm, (40) 

for density n = 0.375 fm“^. One could go beyond this 
limit, speculating that the virial series is asymptotic, in 
the sense that the sum of the first terms gives a rea¬ 
sonable approximation, in spite of the divergence of the 
series, but this is rather risky. 

In Fig. [5]-right, we show, for m = 2,..., 9, the limiting 
values of R for which the accuracy of the truncated series 
is two and ten percent, respectively. For R = 0.3 fm, 
the low density approximation is good enough to get two 
percent accuracy. For R = 0.5 fm, in order to get this 
accuracy, one would need to include two more terms in 
the virial expansion, while for ten percent accuracy, still 
one more term is needed. 

To further improve the equation of state of a gas with 
hard core repulsion, it was proposed , to consider the 
equation of state in the following form 


nT 1 — nv{T) 

This is equivalent to putting 5^ = 1 for all k, instead of 
neglecting them for fc > 2, in the virial expansion (15^ . 
Since 63 = 0.625, the error in the 0{z^) term is reduced 


by almost a factor of two. The errors on the higher or¬ 
der terms are greatly increased, but for sufficiently small 
values of 2; this is unimportant. Comparing with the par¬ 
tial sum Si 2 {z), one finds, that the improved low den¬ 
sity approximation is better than the standard one for 
z < 0.312, i.e. for R < 0.37 fm. 


V. SUMMARY AND CONCLUSIONS 

We have discussed thermodynamics of the excluded 
volume particles at finite temperature and chemical po¬ 
tential, in the low density approximation. Assuming 
Boltzmann statistics, the influence of the excluded vol¬ 
ume on an ideal gas thermodynamics has been derived 
from the change of the free enthalpy and by using the 
Maxwell identities. 

The calculation of the relevant thermodynamic func¬ 
tions and parameters of the excluded volume gas, in 
terms of the two-particle excluded volumes and of the 
thermodynamic functions and parameters of the ideal 
gas, splits into two steps: 

At the first step, it requires calculating the two- 
particle excluded volumes, which is rather standard, 
but one should keep in mind, that the conversion of 
the momentum dependence of the two-particle excluded 
volumes into their temperature dependence, works only 
in the low density approximation and under Boltzmann 
statistics. The single particle exclusive volume is a nat¬ 
ural concept for the single-component gas. For multi- 
component systems, analogous excluded volumes appear, 
but there are ambiguities in their definitions and they 
have no simple geometrical meaning. 

The thermodynamical calculation proposed here, 
which is the second step, is valid in the low density ap¬ 
proximation. In order to go beyond this limit, it would 
be necessary to correct for possible overlaps of excluded 
volumes. 
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The physical picture behind the thermodynamic calcu¬ 
lation, introduced here, is based on the observation, that 
by increasing excluded volume v(T) of particles, the work 
done on the system at constant temperature (T), pres¬ 
sure (p) and particles number (N), equals the increase 
of the free enthalpy. Thus e.g., for the single-component 
gas, the free enthalpy stored in the excluded volumes is 
Npv. Its differential, yields the corrections to the en¬ 
tropy, volume and chemical potential of the gas, and 
through Maxwell identities, provides the relation between 
thermodynamics of the extended and point-like particle 
gases. This new approach, is also valid for the multi- 
component systems. 

A problem of phenomenological importance is that of 
the validity limits for the low density approximation in 
heavy ion collisions. Assuming that all particles are hard 
spheres with the same radius i?, we have shown, based 
on the virial expansion of the equation of state, that for 


R = 0.3 fm the approximation, for the relevant range of 
parameters, is very good. For R = 0.5 fm it introduces a 
significant error, and for R = 1 fm, it is doubtful. 
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